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Abstract 

We revise the conventional nuclear effects of Fermi motion, binding 
and pionic effects in deep inelastic lepton scattering using a relativistic 
formalism for an interacting Fermi sea and the local density approxima- 
tion to translate results from nuclear matter to finite nuclei. In addition 
we also consider effects from p- meson renormalization in the nucleus. 
The use of nucleon Green's functions in terms of their spectral functions 
offers a precise way to account for Fermi motion and binding. On the 
other hand the use of many body Feynman diagrams in a relativistic 
framework allows one to avoid using prescriptions given in the past to 
introduce relativistic corrections in a non relativistic formalism. 

We show that with realistic nucleon spectral functions and meson 
nucleus selfenergies one can get a reasonable description of the EMC 
effect for x > 0.15, outside the shadowing region. 
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1 Introduction 



The EMC effect is probably one of the topics in the interplay of nuclear 
and particle physics which has attracted more attention. Pioneering work 
on Fermi motion and pionic effects [Q, ^, M was followed by many different 
ideas like binding of the nucleons in the nucleus, multiquark cluster effects or 
Q 2 rescaling (see refs. ||, |], [7| for reviews on the topic). Here we shall pay 
attention only to conventional nuclear degrees of freedom, mesons (ir and p) 
and nucleons. One of the interesting ideas along these lines was the effect of 
the nuclear binding || [9|, [TCfl , which, with ups and downs, has come to be 



accepted as, largely or at least partly, responsible for the depletion of R{x) 
2F 2 a{x) / AF 2 d{x) in the region of the minimum. 
Criticism on this latter works was raised in 



11 , where it was shown that 



the introduction of relativistic corrections in the usual nuclear nonrelativistic 
treatment of the binding effects resulted in a flux factor which reduced the 
conventional binding effects [ IT , |12[ . This flux factor led to a different normal- 
ization of the spectral function which preserved the baryonic number calculated 
relativistically fill]. The argument of the normalization of the baryonic number 
is an important one and the idea has met with followers IT! 



T3l[l4||. However, 



it is a prescription on how to convert the nonrelativistic nuclear wave func- 
tion into a relativistic spectral function and the prescription is not shared 
by others ||. Furthermore, these are not the only relativistic corrections as 
shown in [T6|, [T7| . 

This issue justifies a work like the present one, where we construct from the 
beginning a relativistic nucleon spectral function and define everything within 
a field theoretical formalism which uses the nucleon propagators written in 
terms of this spectral function. The relativistic formalism is taken from the 
beginning and the baryonic number is naturally well normalized. Since all the 
nuclear information needed is contained in the nucleon spectral function, one 
does not need to use ordinary nuclear wave functions, which are anyway static 
(no spreading in the energy distribution, which is just concentrated in the 
single particle energies of the shell model) and hence one does not need any 
step to introduce relativistic effects into the nonrelativistic wave functions, as 
done in jy]]. Furthermore, the use of non static spectral functions is important. 
This was already seen in refs. 



13, 



which showed that the use of more 
realistic spectral functions accounting for nuclear correlations resulted in an 
enhancement of the binding effects. The reason is that, for a given average 
binding energy, the approach with a realistic spectral function leads to a larger 
kinetic and potential energy in absolute value than the shell model approach. 

The other issue that we have revisited here is the pionic contribution. Large 
effects from the pionic cloud associated with the pion excess number in the 
nucleus were found in 0, |3|, [?|. We have taken up the idea and recalculated 
these effects within the many body field theoretical approach using input which 
has been checked in a variety of nuclear reactions testing real and virtual pions: 
pionic reactions 0], muon capture [jl9| , inclusive neutrino scattering and 
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photonuclear reactions [21|. In addition we have also included corrections from 
the modification of p-meson cloud in the nucleus, in complete analogy with 
the pionic contribution. 

Altogether we find that the use of the spectral function for the nucleons 
together with the mesonic effects can approximately account for the measured 
EMC effect. 

The calculations are done using the spectral function for nucleons in nuclear 
matter, followed by the local density approximation. This was shown to be an 
excellent tool to deal with photonuclear reactions in the absence of shadowing 
effects |yj . Hence, a natural limit of our results is the region of shadowing, x < 
0.15, where indeed there are discrepancies with the data, and other ingredients 



should be considered that we do not want to tackle [22, 23 



2 Relativistic nucleon propagator in nuclear 
matter 

2.1 Nonrelativistic nucleon propagator and spectral func- 
tions 

Let us recall first the nonrelativistic nucleon propagator for a noninteracting 
Fermi sea, which is given in momentum space by 

G(p°,p)= n 1 '^ . + n(pl 



P° ~ e (Pl + i e P° ~ e (Pl ~ 
where n(p) is the Fermi occupation number n(p) = 1 for |p| < fcp, n(p) = for 
\p\ > hp and s(p) is the nonrelativistic nucleon energy. Eq. (1) can be recast 
as 

G(p°,p) = - Q }— — + 2mn(p) 5(p° - effl (2) 

P ~ £ {P) + %e 

which separates the propagator into the free propagator and the medium cor- 
rection. 

For an interacting Fermi sea the nucleon propagator can be written in 
terms of its nonstatic selfenergy S(p°,p) 

which can be rewritten in terms of the spectral functions for holes and particles 

as 



24 



i* S h (u,p) j , S p (u,p) 



-oo p u — uj — ie Jfj, p u — to + xe 



with the following relationships 
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Sh(p°,p) 



1 JmE(p°,p) 

7r \pO-e(p)-ReT,(p" ,p)] 2 +[ImT, (p° ,p)] 2 

forp < jj, 

_1 JmE(p°p) 

for p° > /i 



let us also recall that the momentum distribution of the nucleon in this inter- 
acting Fermi sea is given by 



ni(p) = / S h (uJ,p) doo 

J — oo 



1 - nj(p) = / S p (uj,p) duo (6) 
with the automatic sum rule 



Sh(uJ,p) duo + / S p (uj,p) duo = 1 (7) 

' — OO w /X 

In passing we also note that in physical reactions nj(p) does not factorize 
out in the physical cross sections because other factors dependent on uo and 
p appear simultaneously in the formulae, and restrictions due to energy and 
momentum conservation do not allow the infinite ranges in the uo integration 
required in eq. (6). Failure to realize that, and the naive substitution of n(p) 
in eq. (1) by rii(p) of eq. (6), as sometimes done, leads to erroneous results 
which can be off by three orders of magnitude in some cases . This gives us 
a warning that we should express all our magnitudes in terms of the spectral 
functions, not the momentum distributions. 



2.2 Relativistic nucleon propagator and spectral func- 
tions. 

The free relativistic nucleon propagator is given by 

p'+M _ M \ ErU r (P)U r (p) Er Vr (~P)^r (~P) \ (g \ 

p 2 — M 2 + ie ~ E(p) \p° — E(p) + ie p° + E(p) - ie J U 

where we have separated in the second member the contribution from positive 
and negative energy states PEfl . M,E(p) in eq. (8) are the nucleon mass and 
the relativistic nucleon energy (p 2 + M 2 ) 1 / 2 , and u r (p), v r (p) are the ordinary 
spinors which we take normalized as u r (p)u r (p) = 1. We recall that u r (p) 
are functions of three momentum and they will be the only spinors which will 
appear in our framework. 

The relativistic nucleon propagator for a noninterating Fermi sea is easily 
derived and, by analogy to eq. (2), can be written as 
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t/) 4- M 

G(p°,p) = 2 T, 2 , • + 2nin(p)(p / + M) 6 (p°) 5 (p 2 - M 2 ) (9) 
which by means of the identity of eq. (8) can be recast as 



G(P°,P) = J^j jl>r(pK(p) 



n{p) n{p) 



p° — E(p) + it p° — E(p) — it 



E r vr(-P)vr(-P)] (1Q) 

p° + E(p) -ie J 1 J 

Apart from the negative energy contribution, which will play no role in our 
problem, the only difference between eq. (10) and the nonrelativistic prop- 
agator of eq. (1) is the presence of the factor M/E(p) and the projector 
over the space of positive energies Y.r u r(p)^r(p), which are both unity in the 
nonrelativistic approximation. 

Now we proceed to construct the relativistic propagator in the interacting 
Fermi sea. We wish to sum the Dyson series for the diagrams shown in fig. 1, 
where although not shown, one would also have other sources of nucleon self- 
energies. We will write them in terms of the operator Y,(p°,p). It will become 
clear later on that we only need the imaginary part of the nucleon propagator 
for the positive energy states, in which case we neglect from the beginning 
the negative energy states (their weight becomes negligible compared to the 
singular part of the positive energy propagator). Hence, for the purpose of the 
present problem, the nucleon propagator needed will be 

M Ur (p)u r (p) ^ M Us (p)u s (p) 
E(p)^p°-E(p) {P ' P> E(p)^p°-E(p) + - 

M U r (p)U r (p) 



E(p) V P° - E(p) - u r (p)^,p)u r (p)^ 

where we have used the fact that £ should be diagonal in spin for spin saturated 
matter which we only consider. 

Comparison of eqs. (11) and (3) shows again the differences between the 
relativistic and nonrelativistic propagators. 

The structure of eq. (11) allows one to define a spectral representation of 
the nucleon propagator by means of Sh(w,p) and S p (u,p) as 



/" 

J— a 



S h (uj,p) 

dui 



p° — uj — irj 
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with the relationships 



S h (p°, P ) 



1 



7T 



+ 



du- 



p° — uj + irj 



'pO-E(p)-^ReJ:(p^p)} 2 + 



for p° < \i 



(12) 



S P (P°,P) 



^ImE(p°,p) 



7T 



for p° > p 



(13) 



where for simplicity E is now -uS-u which is independent of the spin. Eqs. (13) 
are now the generalizations of eqs. (5) using relativistic kinematics. Note that 
S p and Sh defined in eq. (12) are not the nonrelativistic spectral functions 
normally used. Hence one should not expect the same normalization as in 
|TT| , |T2| The normalization of Sh, which we will need, is easily obtained 



by imposing baryon number conservation, as done in |Tl[ . For this purpose we 
evaluate the electromagnetic form factor at q — 0. For the case of the nucleon 
(fig. 2a)) we have (assume all baryons have charge unity for normalization 
purposes) , 



< N\B^\N >= u{p)^u{p) = Bj^; B = 1, p" 



(E(p),p) (14) 



For the case of nucleons in the medium we must evaluate the many body 
diagram of fig. 2b). 



d 4 p 



V iTr[G{p\p)Y}e ip%1 



(15) 



where exp(ip°7]), with rj —>■ + , is the convergence factor for loops appearing 
at the same time and V the volume of our normalization box. 

By means of eq. (12) we can see that the convergence factor limits the 
contribution to the hole spectral function and we get 



< AlB^A >= V J 



d?p M 



(2tt) 3 E(p) 



S h (u,p) du 
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(2vr) 3 E{p) 



(0 + M) on she u % 

2M 7 



S h (u,p) duj 



-"ffam^L^"'^ (16) 

where in the last step we have imposed that this matrix element gives the right 
current with B baryons, in analogy to eq. (14), and P% is the momentum of 
the nucleus. Note that p^ n shell appears in eq. (16) because the operator 
M) onshe u comes from u r (p)u r (p) which depends only on p (it corresponds 
to free particles with p^ = (E(p),p)). Obviously eq. (16) must be evaluated in 
the rest frame of our Fermi sea where all magnitudes are defined. Only p, = 
is then relevant and we obtain the desired normalization 

2V f ^ f S h (u,p)du = B (17) 



(27T)* 

(Note that the factors p °™^» have cancelled in eq. (16)). The factor 2 
is a spin factor. By a simple inspection of eq. (17) we can see that with the 
definition of the spectral functions in eq. (12) one can use eq. (6) to determine 
momentum distributions in both the nonrelativistic and relativistic cases. 

In our formalism we do not have a box of constant density, but elements 
of volume d 3 r with local density p p (r) , p n (r) , the nuclear proton and neutron 
densities at the point r. Hence our spectral functions for protons and neutrons 
are functions of the local Fermi momentum 

k F>p (r) = [tofpptf] 1 ' 8 ; k F>n (r) = [37r 2 P n(r)] 1/3 (18) 
and then the equivalent normalization to eq. (17) is 

d 3 p 



(2< 



3 



S h (u } p } k Fp ,n(r)) du = p p ,n(r) (19) 



In practice we shall work with symmetric nuclear matter of density p(r). Hence 
we have a unique Fermi momentum defined as k F {f) = [37r 2 p(r)/2] 1//3 and then 
one has 



d 3 p 
(2^)3 



S h (w,p, k F {r)) du = p(r) (20) 



— oo 



or equivalently 



/ d 3 r A J 0- £^ S h (u,p, k F (r)) du = A (21) 

with A the mass number of each nucleus. The density p(r) of each nucleus 
is taken from experiment in our case, and expressed in terms of a two Fermi 
parameter distribution for medium and heavy nuclei p7 and mod harmonic 



7 



oscillator for light nuclei |7|, ||. Eq. (21) is fulfilled at the level of 2 - 3 % in 
our case, in spite of the non trivial structure of the spectral function and the 
integrals involved. The small numerical deviation from the right normalization 
is taken care by dividing by the integral of eq. (21) instead of by A in the eval- 
uation of R{x), since a similar integration weighed by the structure functions 
appears in the numerator, as we shall see. Although one could separate the 
contribution of protons and neutrons in the calculation, we have only applied 
the results to nuclei with N = Z, or very close, like 56 Fe and hence, we work 
with the symmetric nuclear matter version. 



3 Deep inelastic electron scattering from nu- 
clei 

Let us recall the basic ideas in deep inelastic scattering. Consider the (e, e') 
process on a nucleon of fig. 3a). The invariant T matrix for the process is 

- %T = ieuJk'WuJk) ^^(-ie) < X\J V \N > (22) 

q 2 

where < X\J V \N > is the invariant matrix element of the hadronic current. 
The cross section for the process eN — > e'X is given in Mandl and Shaw 
normalization |29| by 

1 2m 2M r d 3 k' 2m 

a 



Vrei 2E e (k) 2E(p) J (2vr) 3 2E e (k>) 

6 4 (p + k-k'-Zt lP ^ (23) 

where m is the electron mass, / stands for fermions and b for bosons in the 
final state X. The factor ^Eiv) becomes if we study the cross section on a 
pion. The index i is split in l,j for fermions and bosons respectively. 

In the nucleon rest frame one can then write the differential cross section, 
with Q', E' referring to the outgoing electron, as 

^ 4 TT L> W*> (24) 



dtt'dE' q 4 k ^ 
with a = e 2 /AiT and L' the leptonic tensor 

L'^ = 2kX + K k » + <? V (25) 
and W 1 ^ the hadronic tensor defined as 

W'^ = — W» u (26) 
2ir 
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with 

< X| J^|ff >*< X| J" |ff > (2tt) 4 5 4 (p + g - E^jpJ) (27) 

where q is the momentum of the virtual photon, s p the spin of the nucleon and 
Si the spin of the fermions in X. 



Lorentz covariance and gauge invariance allow one to write W'^ as [|30 



= {*f - r) W X + (^-f [f -^)| (28) 

where W%, W2 are the two structure functions of the nucleon and are functions 
of q 2 ,p.q. 

Now we evaluate the cross section for (e, e') on the nucleus. In order not 
to miss flux factors and be able to write everything in terms of propagators, 
we evaluate the electron selfenergy corresponding to the diagram in fig. 3 b). 
We obtain 

Z^R n ^(g) Z^fi (29) 
which for unpolarized electrons can be written as 

with IP"(g) the photon selfenergy. 

In order to evaluate the cross section we need only ImL{k) which can be 
easily evaluated by means of Cutkosky rules [26| . 



E(Jfe) -> 2iImZ(k) 

D(k') — > 2i6(k'°) ImD(k') (boson propagator) 

n^(g) -> 2«0(g°) Jmn^(g) 

G(p) — > 2i8(p°) ImG{p) (fermion propagator) 



(31) 



with the result 



ImYiCk) = e 2 [ 1 %°) ImW u (q) - A — U (32) 

1 J J (2tt) 3 2E e (k - q) q A 2m ^ 1 J 

with q° = k° - E e (k-q). 
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The cross section is readily evaluated from there. Inspection of eq. (11) 
for the relativistic fermion propagator tells us that the electron width is given 
by 

2?T7 

r(jfc) = ImX{k) (33) 

E e (k) 

by means of which the contribution to the cross section from an element of 
volume d 3 r in the rest frame of the nucleus is 

da = TdtdS = T— dldS =-d 3 r = 
dl v 

„ EJk) ,o 2m o . , 

= T d 3 r = — — Im £ d 3 r (34) 

Hence we immediately write the (e, e') cross section in the nucleus as 

with q° = k° — E e {k — q) = k° — k'°, which is always positive in this experiment, 
so we drop the 9 {(f) function. Comparison of eq. (35) with eq. (24) used for 
nuclear targets tells as that 

W'F(q) = -4 - Id 3 rim W v {q) (36) 

Next we evaluate IP" {q) corresponding to the right hand side of the dia- 
gram of fig. 3b) using again the Feynman rules in terms of propagators. We 
have 

d4 P ^ ^ ^ ttTV f d4 Pi 



« = (-)/(2^ iG W E - E - n -/(2.)< 



UiiGMUj iDj^X-ife 2 < X\J»\H >< X\J V \H >* (2tt) 4 S^q+p-Eg^ 

(37) 

which by means of Cutkosky rules (31) and the use of free propagators, eq. 
(8), for the final states and the medium propagator, eq. (12), for G(p), plus 
eq. (36), leads immediately to 

i- s x s Sp s Si nf =1 J ^Uu f (j^j u j£b (JL) 

< X\J»\H >< X\r\H >* {27r) 4 5 4 {q+p-^ = iP'i) (38) 
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which by means of eq. (27) can be rewritten for an isospin symmetric nucleus 

as 



with 



p = (p°,f); W'jT = \{W'r + WT) (39) 

Note that W'^ip, q) appears with the off shell arguments of p, the bound 
nucleon. 

In the steps from eq.(37) to (38) the spinors u(p) are included in the matrix 
elements of the currents and we have considered that there is necessarily a 
fermion loop (hence the first minus sign in eq. (37)) with a free particle 
in the final state and the nucleon in the medium in the initial state. The 
corresponding energy integration in the loop (if we had used a Wick rotation 
explicitly instead of Cutkosky rules) necessarily picks up the hole part of the 
propagator of eq. (12), and at the same time relaxes the condition 9(p°) of 
Cutkosky rules which does not appear for the hole part. 

In eq. (39) there is an apparent lack of normalization, since assuming 
W'^ u constant (which actually cannot be in practice) we would expect Wjf" = 
AW ,flu . However given the normalization of the spectral function in eq. (21), 
this is not the case. In eq. (39) we get the extra factor which does not 
appear in eq. (21). It is interesting to see the meaning of this factor in eq. 
(39). If we look at the formula of the eN cross section in eq. (23), and by 
means of eqs. (25) and (27) we find 

a 2 M r d 3 k' 1 



v rel E e (k)E(p) J E e (k') q 4 

where L' W'^ is a Lorentz invariant and the content of the integral in eq. 
(40) also. For collinear frames of reference we also have 

v rel E e {k)E{p) = Mk (41) 

where k is the electron momentum in the frame where the nucleon is at rest 
and hence 



i4kh L ^ wT (42) 



°~ k J E e (k') q 4 

Now if we have a system of moving nucleons the cross section for scattering 
of the electron with the nucleus can not be obtained as a sum of individual cross 
sections, because the relative eN flux is different for each nucleon. Instead, 
one has to sum the probabilities of collision per unit time for each nucleon and 
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divide by a unique flux, the one relative to the CM of the nucleus. By taking 
for v re i in eq. (40) the velocity of the electron with respect to the CM of the 
nucleus and suming over all the nucleons in eq. (40) we will be calculating 
the electron nucleus cross section. Hence we obtain in the rest frame of the 
nucleus 

- Tl-Fm7 L '" 4 J" 3r J (ft m L w, p) wtm 

(43) 

where k is the electron momentum in the nucleus rest frame. Since our nuclear 
cross section is given in terms of Wjf v by (see eqs. (35), (36)) 

°? f jj w ,»v (44) 

k J E e (k') ^ A ^ ' 

then eq. (39) follows immediately. 

The previous discussion has shown that the factor -g^L is a factor appearing 
in the probability of reaction per unit time for each nucleon, and remains in 
the integral when we divide by a unique flux in order to obtain the nuclear 
cross section. It is thus a Lorentz contraction factor. 

In the limit of small densities, when M/E(p) = 1, eq. (39) with the 
consideration of eq. (21) would give Wj£ v = AW'^ V as it should be. Eq. (39) 
accounts for Fermi motion and binding and includes the relativistic Lorentz 
contraction factor M/E(p) and the change of the arguments in W ltiV {p, q). 
Note in passing that the relativistic factor m/E e (k) of eq. (33), which we 
extracted from the relativistic propagator of eq. (11), has been essential to 
provide the right normalization. 

4 Contribution from the pion cloud 

Let us first see the free pion structure function. The same formula eq. (24) 
is used for pions and this defines W^ v . Given the normalization of the fields 
which we follow |[29| , the cross section of eq. (23) contains the factor l/2u(p) 
instead of 2M/2E(p). Hence, this means that the definition W!^ u is given by 
eq. (27) dividing the right hand side of the equation by 2m,,- (and obviously 
the average over the spin of the nucleon does not appear now for the pion 
case) . 

In order to derive the contribution from the virtual pions in the medium 
we evaluate again the electron selfenergy related to the diagram of fig. 4. We 
can save all the steps given before simply by noting the differences in the two 
cases: 
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i) The bound nucleon propagator is substituted by a pion propagator. From 
the use of Cutkosky rules we must change 




— oo 



duS h (u,p) 5 (p 



,o 



by 



26(p")ImD{j>) 



(45) 



with D(p) the pion propagator (in the medium). 

ii) One must take into account that W!^ u is given by eq. (27) divided by 
2m„-. 

iii) The sum over spins of the bound nucleon in eq. (37) does not appear 
now for the case of the pion. 

iv) There are three charged states of pions. 
With only these four rules we can already write 



Now there are two obvious subtractions to eq. (46). First one should 
subtract the contribution from a free pion, which has nothing to do with 
medium effects. However, this is zero because one electron can not decay into 
another electron, one pion and X. Then assuming the pion is dressed in the 
medium by exciting ph and Ah as we shall do, one is left with the contributions 
shown in fig. 5. 

Now there is no problem to get a contribution since ImD(p) gets strength 
fromp/i excitation. The physical channel would correspond to e — > e' + X +ph, 
or equivalently eiV — > e'N'X, which is now allowed. The physical channels are 
easily visualized by cutting the intermediate states in the diagramas with a 
horizontal line and placing on shell the particles cut by the line. This is actually 
the essence of Cutkosky rules to obtain the imaginary part of the selfenergy 
of a diagram. 

The former discussion also tells us that part of what we are calculating is 
already contained in the nucleon structure function. This is because we are 
also calculating the contribution from the pions contained in a free nucleon. 
This has to be subtracted. This is easily done by substituting in eq. (46) 



since we substract A times the contribution from the pion cloud to the struc- 
ture function of the free nucleon. In technical words, we can say that we are 
only considering terms with at least two ph or lphlAh (the most important 
terms) etc. in fig. 5 (up to Pauli blocking corrections in the ph Lindhard 
function, which are automatically included by the procedure of eq. (47)). 
Hence the genuine pionic contribution is given by 



W'Z = 3 / d 3 r 




ImD{p) — > 5ImD{p) = ImD{p) — p 



dImD(p) 
dp 



p=0 



(47) 
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W'X: = 3 j d 3 r J-0- 6(p°) (-2)5ImD(p) 2m. (p, q) (48) 

Eqn. (39) and (48) are the basic equations which provide the nucleonic and 
pionic contributions. 

In passing we can mention that the distribution of the excess number of 
pions, per unit volume in the nucleus, 5N n (p), which contains the averages of 
< a+ a p >, < cip al p > and < a p a_ p >, is given by 

snap! _ _ 3 r f_ SImD{p) m 



2uj{p) Jo 2n 

such that in the case of a structure function W^ u (p, q) independent of p°, eq. 
(48) could be considered as a weighed integral of the pion structure function 
over the pion excess distribution in the nucleus. However, the strong depen- 
dence of W^ u (p,q) on p° (imposed by energy and momentum conservation) 
does not allow that factorization, and hence a relationship of the pion ex- 
cess number with the pionic contribution to the structure function cannot be 
established. The apparent extra factor 2 which we obtain in this counting 
(apart from the Lorentz contraction factor, m- K /uj{p)) is obtained because one 
is automatically accounting for the imaginary part of the Compton 77? am- 
plitude which is crossing symmetric and contains the two diagrams of fig. 6, 
while the pion structure function for on shell pions contains only the imaginary 



part of the diagram 6a). (see refs. [31] and [B3] for an elaborate discussion 



of these issues in the problem of the pion cloud contribution to K + nucleus 
scattering). It is worth noting that there is no overlap between the Feynman 
diagrams accounted for in the pionic contribution, fig 5, and those of the nucle- 
onic contribution with Fermi motion and binding, which come from selfenergy 
insertions in the nucleon line of fig. 3b). Hence, these contributions to the 
nuclear structure function are independent. 

5 The Bjorken limit 

We have evaluated the contribution of nucleons and pions to the hadron struc- 
ture function of the nucleus. We now proceed to write these expressions in 
terms of the Bjorken structure functions [|33|| . 

For nucleons (and similarly for pions or the nucleus) one introduces the 
Bjorken variables 

X= 2p^^=M ]Q= - q (50) 

and for large values of q° and Q 2 simultaneously and fixed x one has the 
Bjorken scaling 

vW 2 (x } Q 2 ) =F 2 (x) 
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MW x (x,Q 2 ) = ^ (a;) (51) 
and the Callan-Gross relation 

2xF 1 (x) = F 2 (x) (52) 

up to some, QCD corrections in InQ 2 . Since the same Q 2 will be chosen for 
the nucleus and the nucleon and we perfom ratios of structure functions we 
shall not worry about this dependence here. 

In view of the relations (51), (52) the most practical way to proceed is to 
work with transverse components of W ,fiu . For this purpose assume q along 
the z direction, as usually done in the study of the e, e' reaction, and evaluate 
W' xx . We find from eq. (28) 

1 M 2 2 ~ M M 2 v 



F\(x) 

= in the Bjorken limit (53) 

this component has the virtue that the coefficient of W± is independent of p 
and hence is the same for on shell or off shell nucleons, or pions, or the nucleus. 
Hence we can write 



FiaAxa) -4 [d*r [*r 4- dp°s h (A P ) Fin{xn) 



M A J J (2tt) 3 E(p) J-oo * ^ M 

FiaA x a) _ o f , 3 f d4 



Sjd 3 r J 0(p°)(-2) 5ImD(p) 2m n 



M A J J (2tt) 4 * " > v m ^ 

x A = -q 2 /2M A q° = - , with x = -q 2 /2Mq° 

xn = -q 2 /2pq ; = -q 2 / - 2pq (54) 

where the extra minus sign in x n is because of the direction of p in fig. 4. 

Here F 1N (x) = (F lp (x) + F ln (x))/2 as implicit in eq. (39). 

We still have to exert some caution since < xn < 1- On the other hand 
x n < 1 but x-x > x because in our scheme the emerging particle from the 
coupling of the pion to a nucleon of the nucleus is on shell when we excite 
ph, Ah with the pion and take the imaginary part of Dijp) 0. 

Since usually one compares ratios of the F 2 structure functions, this is 
easily accomplished by making use of the Callan-Gross relation (52) and we 
find 
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d * r / To - Tv^ / dp S h ( P °, P ) — F 2N (x N ) 9{x N ) 9(l-x N ) 

J {Ali) A F(p) J-oo Xj^ 



(55) 



F 2 aA x a) = -6 J d 3 r J -0^ 9(p ) 5ImD(p) ^ 2M F 27T (x^) 9{x lx -x) 9(l-x*) 

(56) 

where we have again F 2N = (F 2p + F 2n )/2 as implicit in eq. (39). 

In the Bjorken limit the evaluation of eq. (55) does not require the knowl- 
edge of the variable q since 

x p° — p 3 x —p°+p 3 

x~^^ M ' x^~ M ^ ' 

but it is implicit in the structure functions which are taken at a certain Q 2 . 

When reaching this point it is worth considering also the contribution of 
the p meson cloud. Both the pion and the p meson couple to nucleons and 
delta with derivative couplings, which give rise to relatively large meson self- 
energies in the range of momenta which contributes to the structure functions. 



Furthermore, as found already in ref. [31|, the pion cloud contribution comes 



mostly from the combined ph and Ah excitation, with the ph on shell (in 
SImD(p)), and one finds a negligible contribution of two ph excitations. Only 
it and p can excite the Ah components and this makes these two mesons special 
when looking at the mesonic contribution to the nuclear structure function. 
In addition there is an interplay between 7r and p exchange. Indeed, the p 
meson, through nuclear correlations, contributes both to the longitudinal and 
transverse parts of the spin-isospin ph and Ah interaction and it is an impor- 
tant element contributing to the Landau-Migdal g' parameter in a microscopic 
derivation of this interaction |34| . The value of this parameter governs to some 
extend the pionic (and p- meson) contribution to the structure function [[|. 

The contribution of the p-meson cloud to the structure function is given, 
by analogy to eq. (56) by 

F 2A , P M = -12 / d 3 r [ ^-9(p°)5ImD p (p)-2MF 2p (x p ) 
J J {^) x p 

9(x p - x)9{l - x p ) (58) 
where D p {p) is now the p-meson propagator and F 2p (x p )is the p-meson struc- 



ture function, which we take equal to the one of the pion following refs. |35]. [36 
In addition x p is also given, in analogy to eq. (57), by 

- = < 59 > 

x p M 
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Eq. (58) contains an extra factor of two compared to the pionic contri- 
bution of eq. (56). This is because of the two transverse polarizations of the 
p-meson and the fact that the coupling of the p to nucleons and deltas which 
we consider, following ref. |34| is only of transverse nature ((<? x p)e for nu- 
cleons and (S + x p)e for deltas, with e the polarization vector of the p meson 
and S + the spin transition operator from spin 1/2 to 3/2). 

The expression of Fia,n{ x a) in eq. (54) shows the same lack of normaliza- 
tion discussed in connection with eq. (39), since assuming F\n(xn) constant 
(which is not the case) Fia,n is not A times Fin, due to the Lorentz contrac- 
tion factor -g^. The same can be said about eq. (55) in general. However, 
eq. (55), shows a particular normalization property. Indeed, if we take an 
ensemble of nucleons on shell (p° = E(p)) and x = (and hence x^ = 0, by 
virtue of eq. (57)), then the Lorentz contraction factor -^Mry cancells in average 
the dynamical factor and we get F 2 a,n(0) — AF 2 n{0) for on shell nucleons. 
This cancellation, however, will not show up when the nucleons are off shell 
since in our formalism we still obtain the factor but the dynamical factor 
— will now be different. 

X N 

In this discussion we are implicitly assuming that eqs. (54) and (55) stand 
as they are for the case of off shell nucleons simply by taking for x^ the 
expression of eq. (57) using the nucleon off shell variables. This is certainly 
the easiest form of the analytical continuation in the off shell regime, although 



there have been other prescriptions in the Literature [[37], 0, [39[ . 

We would like to justify our assumption. Indeed, in our frawework we use 
nucleon propagators which are based on the free spinors (eq. (12)). Hence in 
eq. (38) the matrix elements would be defined in terms of free spinors, the 
final particles are free particles and the off shell dependence appears only in 
the S( ) function. Obviously one can not look in detail at all channels implicit 
in eq. (38). However, one can use the same philosophy in the parton model 
which is used to find out the scaling of the structure functions. We follow 



here the steps of ref. |33| and assume that partons carry a fraction x q of the 
nucleon momentum p and its mass M, that the electron parton amplitude is 
given by the on shell expression and that the outgoing parton is a free one, 
but the S( ) function appears with the off shell variables. 
One finds then 

W x = W' xx = J2j dx q f t (x q )e^(-q 2 )-^6(q 2 + 2pqx q ) 

= Ee*/<(^)47 = ^T^W = — (60) 
Y 2M M 2pq K J 

where F\ appears whith argument x^ which is defined in terms of the off shell 
variables. Similarly one obtains 

V A W2 = Y / e 2 f i (x N )x N = F 2 (x N ) (61) 
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whith p.q, and denned in terms of the off shell variables. 

We would like to note here that our finding F 2 a(0) ^ AF 2 jv(0), or equiv- 
alently R(0) ^ 1, where R(x) = F 2 a(x)/AF 2N (x), is not so unconventional. 
Indeed, in ref. |16| the authors also find this property and they literally quote 
"Note that a feature of the results is that R(0) 7^ 1... It does not reflect 
any violation of baryon conservation, which is ensured by the normalization 
condition", in our case eq. (21). 

The issue of the normalization still stirs much controversy. We have de- 
voted many thoughts to it throughout this paper looking at it from some 
points of view not discussed before. We see that for an ensemble of uncorre- 
lated nucleons F 2 a,n(0) = AF 2 n(0) but as soon as interactions are accounted 
for this normalization is lost. A similar thing would happen should we eval- 
uate the structure function F 3 a, which appears in neutrino scattering. This 
structure function when integrated over x is normalized to 3 for the nucleon if 
one assumes models with only three valence quarks, or to 3A for the nucleus. 
Once again this normalization would be lost if QCD interaction corrections 
are accounted for ||D| and equivalently if NN interactions are included in the 



nuclear case, and experimentally this is the case [40 . Also, experimentally 
F 2 a,n(Q) 7^ AF 2N (0). Clear as the question looks to us, we are aware that this 
point of view is not universally accepted. Further thoughts and discussions 
on the issue should be welcome which would help settle the question in a way 
acceptable to all. 

Eqs. (55, 56, 58) are the final equations which we use in the analysis. 
For the nucleon and pion structure functions F 2N (x), F 27r (x) we take the 



experimental values of refs. [41, 42 



6 The meson propagators 

The pion propagator in the medium is given by 

D(p) = [p 02 -f-ml- U n (p°,p)}- 1 (62) 

with H,,- the pion selfenergy. We consider the contribution of the ph and 
Ah excitations to the pion selfenergy in connection with the Landau-Migdal 
correction in the pionic channel, as well as off shell pion nucleon form factors. 
Since for 5ImD we need D — Dq, it is practical to write 

J^F 2 (p)p 2 II* (p) 
D{p) - D (p) = D 2 (p) / 2 g I ; : ' (63) 

where F(p) is the nNN form factor, which we take of the monopole type 

A 2 — 771^2 , , 

F (P) = ( 64 ) 



18 



with A = 1300 MeV [H and / 2 /4vr = 0.08. V[{p) is the longitudinal part of 
the spin-isospin interaction and II* (p) is the irreducible pion selfenergy, which 
contains all selfenergy diagrams which are not connected by V[{p). 
For the p-meson we can write 

-^C p F 2 (p)p 2 U*(p) 

^-^-^ ^U (65) 



with Vr(p) the transverse part of the spin-isospin interaction, C p = 3.94 [4~3"| 
and F p (p) the pNN form factor given by 



A 2 — rn 2 



and A p = UOOMeV 14 



For V[(p),Vt(p) we take the expressions which are derived from a model 



with 7r and p exchange in the presence of short range nuclear correlations [ji4 
which are given by 

V£(p) =p 2 D (p)F 2 (p) -p 2 D (p)F 2 (p) 



~ \q 2 Mp)F\p) ~ \qlDop{p)F 2 p {p)Cp (67) 



V* T (p) = p 2 D 0p (p)F 2 (p)Cp - \q 2 Mp)F 2 {p) 

-{p 2 + \ql)Do P {p)F 2 p{p)Cp (68) 

Here q c ~ 780MeV is the inverse of a typical correlation distance and 
D(p), F(p), D p (p), Fp{p) are the corresponding propagators and form factors 
substituting p 2 by p 2 + q 2 . The irreducible selfenergy II* (p) is in our case the 
sum of the Lindhard functions Un{p),U 1 \(p) for ph and Ah excitation with 



the normalization and analytical expressions of the appendix of ref. |45[ . It 
is interesting to note that for the values of p°,p which contribute most to the 
structure function, both V[(p) and V^(p) are negative and this leads to positive 
values of F 2 a,-k and F 2 a, p - 



It was also found in []3T| that in order to evaluate the contribution of the 
pion cloud by using the pion propagator, as done here, it is important that it 
satisfies the sum rule 

roo (^rp 

- — 2p° ImD(p°,p) = 1 (69) 

JO TT 

which expresses the equal time commutation relation of the pion fields. Our 
model satisfies this equation at the level of one per thousand. 
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7 The nucleon spectral function 



Section 2.2 has established the framework for the relativistic nucleon propaga- 
tor which we need here. The only input needed is the nucleon selfenergy. We 
take it from the work of ref. K5I . This is a semiphenomenological, quite suc- 



cesful approach, which uses as input the NN cross section and the spin-isospin 
effective interaction. This allows one to evaluate JmE, which is remarkably 
close to ImS of the elaborate many body calculations of ref. |^7|, |48[ . The 
real part is evaluated by means of a dispersion relation, and the Fock term 
from the pionic contribution is also included. Only pieces of the Hartree type, 
which should be independent of the momentum, are missing for which one 
needs more information. Hence, up to an unknown momentum independent 
piece in the selfenergy the rest of the nucleon properties in the medium can 
be calculated, like effective masses, spectral functions, etc, which are also in 
good agreement with sofisticated many body calculations |49], ^0f| . Actually, 



what might appear a drawback is now a welcome feature because since the 
proper binding energy is an important ingredient in the EMC effect, we also 
include phenomenologically a function C(p) in the nucleon selfenergy and de- 
mand that the binding energy per nucleon be the experimental one for each 
nucleus. Then the model is complete, realistic and technically much simpler 
to handle than the sofisticated many body calculations [|9|, [5(| . 

With this improvement, momentum distributions and average binding en- 



ergies are also in good agreement with other infinite nuclear matter ETJ and 



finite nuclei calculations [52] 



A small inconvenience appears because the selfenergy of ref. EB[ is evalu 



ated non relativistically. However we have checked that a proper calculation 
including relativistic factors of the type M/E(p) in the nucleon propagators 
in the integrals over the loops which appear in the evaluation of IrriE in |[46| , 
would only introduce corrections in JmS below the level of 10%. Second, we 
have changed ImH by 10% and found that the ratio R(x) changes only at the 
level of 1%. Thus we take the values for £ from ref. f46[ and use them in 
the relativistic propagators of section 2.2. For the average kinetic and total 
nucleon energy we have 

<T>= jJ d3r Jwf> (E{f> ~ M) I" x Ship °' p)dp ° 
< E >-lI d3r lw? /-> (p °- p)p0 *° (70) 



and the binding energy per nucleon is then given by the sum rule [53 



\Ea\ = - l -(<E-M> < T > ) (71) 



which is also used in p% pq] in connection with the study of the EMC effect. 
We take experimental numbers for each nucleus for \E^\ and adjust the func- 



20 



tion C(p) to fit \Ea\- We take C(p) linear in the density, Cp(r). This quantity 
provides around 30MeV repulsion at p = p in most of the nuclei. Detailed 
values for < T >, < E > and Ea can be seen in Table I. 



As noted in refs. |p~3| , |54| , [55| , the use of nucleon propagators in terms of non 
static spectral functions leads to bigger values of the average kinetic energy 
and | < E — M > | than the shell model of the nucleus and as a consequence to 
reduced values of R(x) (for x < 0.7). We can see this here also by taking the 
uncorrelated Fermi sea and adding a function Dp(r) to the ordinary Thomas- 
Fermi potential Vtf(^) = — kp{r) 2 /2M such as to get the same binding energy 
via eq. (71). We can make use of the same formalism by simply considering 
that 

S^ s (p°,p)=n(p)5(p°-E(p)-J:) 

H{r)=V TF {r)+Dp{r) (72) 

Eq. (72) associates one energy to a given momentum, (the essence of the 
shell model in infinite nuclear matter) while the spectral function has a peak 
around the quasiparticle energy and then spreads out at larger values of p for 
a given energy. This results in a larger value of the average kinetic energy. 

In Table I we show the results for different nuclei and compare them with 



those of ref. |54 ]. As one can see, the results that we obtain with the uncor- 
related Fermi sea and the spectral functions are remarkably close respectively 
to those of the Hartree Fock and spectral function used in []54| . 

We shall evaluate R(x) using both the spectral function approach and the 
uncorrelated Fermi sea. We have integrated over the momentum up to four 
times the Fermi momentum for each energy. This gives the normalization of 
A at the level of 2 — 3% which is sufficient for our purposes, but can lead to 
higher uncertainties in the kinetic energy, which weighs the integral with a 
higher power of p. Even a conservative error of 20% in the kinetic energy has 
repercussions in the EMC effect only at the level of 2%. 

8 Results and discussion 
8.1 Nucleonic contribution 

In fig. 7 we show the results for R(x) from the nucleonic contribution calcu- 
lated as 

Rn{x) = 73 
AF 2N (x) 

We do not divide by F2d{x) as experimentally done. One reason for it is 
that the techniques used here with the local density approximation cannot be 
used for deuterium and hence we cannot calculate the nuclei and deuterium 
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with the same model. However the price we pay is small. We rely upon the 



results for F 2D (x) calculated for deuterium in |Tj|. We can see there that for 
< x < 0.7 F 2D (x)/ F 2N (x) ranges between 0.98 — 1. Hence, our results 
should be increased by 1 — 2% and we will not worry about this amount. More 
serious is the region for x > 0.7, where due to Fermi motion the former ratio 
increases rapidly. Hence, in that region we should expect to overcount the 
experiment as it is indeed the case. 

In fig. 7 we plot the results obtained for 56 Fe. This nucleus has N ^ Z but 
by a little amount. Furthermore, the experimental results are corrected by the 
isoscalarity factor to convert them into an equivalent isoscalar nucleus |U| , 



hence our calculations done for symmetric nuclear matter are appropriate. We 
see a minimum around x = 0.6 as in the experiment and a steep rise at x > 0.75 
as it correponds to Fermi motion [j5"7| , 0, B|. The region around x = 0.6 agrees 
well with experiment. However at x ~ 0.15 — 0.2 we are below the data by 
about 10 — 15%. This region will be filled up latter by the mesonic contribution. 

In the same figure we show the results obtained with the uncorrelated Fermi 
sea (local step function distribution) of eq. (72). We observe that Rn(x) takes 
values closer to unity than the results with the spectral function for x < 0.6 
The reduction of Rn{x) with the use of the spectral function with respect to a 
static picture of the nucleus, like Hartree Fock or the equivalent uncorrelated 
Fermi sea in our case, was already shown and explained in ref. ||54j| . The 



results here are qualitatively similar to those in |54| and the explanation lies 
in the increased binding provided by the spectral functions. Both the factor 
x/xn of eq. (57), as well as the restrictions of phase space 9(xn)9(1 — x^), 
are responsible for the decrease of Rn{x) in this region. 

On the other hand there is a novelty in these results with respect to those 



in [|54[ . Rn(x) does not go to 1 at x = as in [54j] and are systematically 
lower in all the range of x. The reduction at x = is easy to see from the 
formulae. In eq. (55) F 2 n{xn — 0) will take a constant value, and with respect 
to the normalization integral of eq. (21) the novelties in eq. (55) are the extra 
factors M/E(p) and x/xn which both go into reducing the contribution of the 
integrand for off shell nucleons. 

The overshooting of the results in the region of x ~ 0.8 was already an- 
nounced as a result of dividing F 2 a by F 2N and not F 2D /2. But the qualitative 
features due to Fermi motion are reproduced. A very detailed discussion of 
these effects is given in ref. |J, but qualitatively we can see that 

x N M 



x N 



x p° — p 3 



(x -> 1) (74) 

and one can pick up values of p 3 in the integrations such that x^ < 1. Hence 
F 2 a,n will be different from zero while F 2 ^(x = 1) = 0, and Rn(x) from eq. 
(73) necessarily goes to oo. 
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It is interesting to call on the attention to the crossing of the two lines in 
fig. 7. The binding effects reduce Rn(x). On the other hand Fermi motion 
increases Rn(x) close to 1, as we noted. Fermi motion is more important in 
the interacting Fermi sea because now one has larger momentum components. 
On the other hand the interacting Fermi sea has also more binding. As a 
consequence we see that at x < 0.6, where the binding effects dominate over 
the Fermi motion, Rn{x) for the interacting Fermi sea is smaller than with the 
non interacting Fermi sea, while for x > 0.7, where the Fermi motion effect 
dominates, the situation is just opposite. 



8.2 Mesonic contributions 

In figs. 8, 9, 10, 11 we show the mesonic contribution to R(x), together with 
the nucleonic one discussed above, for different nuclei, 6 Li, 12 C, A0 Ca and 
56 Fe. The general features are the same in all nuclei, but, of course, the 
mesonic contribution is smaller in lighter nuclei. The mesonic contribution is 
calculated with the structure function of ref. fH2l. We have also calculated it 



with the older structure functions of ref. 58, 59 . We find that around x = 0.2 



the results for R(x) decrease in about 3% if one uses the pion structure function 



of ref. |58| and increase in about 5% if one uses the one of ref. |5"9" |. This should 
give us an idea of the uncertainties of this contribution. In fig. 11, for 56 Fe, 
we split the mesonic contribution into the pion and p-meson ones. We observe, 
that although the pionic contribution is bigger, the one from the p-meson cloud 
is also important, and both of them are positive in all the range of x. Similarly, 
as obtained in other calculations |2|, |3|, [7j the mesonic contribution vanishes 
around x = 0.6 and increases as x decreases. We can see that thanks to the 
mesonic contribution the agreement with the data becomes much better. The 
slope of R(x) from x = 0.15 to x = 0.6 is not reproduced by the contribution 



of the nucleons alone, a feature which is shared by the results of ref. [Q. The 
mesonic contribution comes to produce the right slope. 

As we said, the qualitative features of the pionic contribution are similar 
to those in 0, §, 0]. There are also some differences. In ref. the pionic 
contribution is evaluated assuming different amounts of pion excess in the 
nucleus, but no evaluation of this excess is made. Furthermore one should 
recall our warnings in section 4 not to use the excess number in the evaluation. 
In ref. [[| an actual evaluation is done of the pionic contribution. Even if 
the formalisms here and there might look quite different, they are actually 
quite similar and one can see that ImD(q) = \D(q) 2 \ ImU n of eq. (48) 
appears in ref. || as Jmri 7r |_D (g)| 2 , with \D(q)\ 2 changed to |-D (g)| 2 , which 
is not very problematic when the pions are off shell. Furthermore, the Ah 
contribution is obtained from an extrapolation of the pionic atom data, which 
is more problematic when one goes to the off shell situations which one finds 
here. Furthermore, this Ah is taken real in ||, while here the Ah Lindhard 
function is explicitly evaluated as a function of q°, q by keeping the A width. 
An accurate evaluation of the imaginary parts of the diagrams is necessary 
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in order to fulfill the sum rule of eq. (69). On the other hand what one 
evaluates in ref. [Q] is the fractional increase of R(x) with respect to the one 
in the Sullivan process (deep inelastic scattering with the pion cloud of a free 
nucleon) ^T| and not the absolute value of the contribution of the pionic 
cloud. 

In ref. J7J the absolute contribution of the pion cloud is obtained by sub- 
tracting the " free " parts contained in the response function of a free nucleon, 
as done here. The formalism is similar to ours but some approximations are 
done which are improved here. For instance the Ah contribution is again taken 
real and other approximations are done to relate some magnitudes to the pion 
excess number (recall our warnings about this). 

The input which we have used for the meson nucleus interaction, Vl(p), Vr(p), 
etc, is the one of section 6. It has the virtue of having been tested in a large 
variety of reactions and we have not changed it here. This gives us much 
confidence about the strength of the mesonic contribution obtained here. The 
fact that it fills up the part of R(x) missed by the nucleonic contribution is 
certainly a welcome feature which reinforces our confidence on this mesonic 
model. 

The agreement of the results that we obtain with the experimental data 
can be considered rather good by comparison with results obtained with other 
theoretical approaches. The trend of the data is well reproduced and the 
remaining discrepancies are not incompatible with the intrinsic theoretical 
uncertainties of our model, particularly of the mesonic contribution and more 
specially of the p-meson cloud which is somewhat tied to the form factors and 
the nuclear correlation function, the fact that we divide by F 2 n{x) instead of 
F 2D (x)/2 ,etc. 

Our results for R(x) in m Ca and 56 Fe look already very similar and we 
have checked that R(x) does not change much for heavier nuclei. Obviously 
one should take into account that for heavier nuclei N ^ Z and our approach 
with symmetric nuclear matter should be less accurate. Actual calculations 
keeping different neutron and proton densities lead to a slight decrease of the 



minimum of R(x) as Z increases [55 



9 Conclusions 

We have evaluated the nucleonic and mesonic contributions to the ratio R(x) 
of the EMC effect, with particular emphasis on an accurate treatment of 
effects shown in the past to be relevant, like binding effects, Fermi motion 
and a dynamical (non static) treatment of the nucleons and the mesons in the 
medium. 

In order to avoid having to take some prescription on how to include rela- 
tivistic effects in the approach, which has led to many discussions in the past, 
we started with a relativistic approach from the beginning and have made a co- 
variant treatment which allows us to write all magnitudes in terms of nucleon 
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and meson propagators in the medium. The approach was made practical by 
evaluating the relevant magnitudes in an infinite medium and calculating the 
structure functions in finite nuclei by means of the local density approxima- 
tion. This procedure is fine for x > 0.1 but certainly breaks down for x < 0.1 
where there is nuclear shadowing. 

We could see that the use of the spectral functions to construct the nucleon 
selfenergy was relevant in reducing somewhat the ratio R(x) with respect to a 
static picture of the nucleus, like a shell model, or in our case an uncorrelated 
Fermi sea. This reconfirmed qualitatively earlier findings in the same direction. 

Although our results for the nucleonic contribution differ somewhat from 
other results in the literature, we share their conclusions that the nucleonic 
contribution alone does not explain the data, particularly the slope from x = 
to x = 0.6. 

On the other hand, we evaluated the contribution from the pion and p- 
meson clouds rather accurately. Recent work done before on the contribution 
of the pion cloud to the K + nucleus selfenergy had taught us some important 
lessons, particularly the importance of using an input which satisfies a sum 
rule, not trivial to satisfy unless the analytical properties of the pion selfenergy 
are strictly fulfilled, and the need to avoid any relationship to the "pion excess 
number". In addition, experience gained in dealing with reactions which in- 
volve real and virtual pions allowed us to use information on the pion nucleus 
selfenergy which is realistic enough and has been tested in many such reac- 
tions. Hence, we consider the present calculation of the pionic effects as an 
improvement over work done in the past and we think these results are rather 
reliable. However, there are still small uncertainties in the pionic contribution 
stemming from different results for the pion structure function obtained in 
different analyses of the Drell-Yan process. 

The strength of the pionic effects is moderate. So is the one from the 
p-meson cloud, but, when they are added to the nucleonic contribution one 
obtains a good description of the data in the region of 0.1 < x < 0.7. 

In summary we could conclude that the main features of the EMC effect 
can be described in terms of conventional degrees of freedom, nucleons and 
mesons. It does not exclude explanations in terms of more elementary degrees 
of freedom like quarks and gluons. It is simply a question of which degrees 
of freedom are more economical and transparent, as stressed by Jaffe in ref. 
15] . The fact that we could deal with these degrees of freedom with a certain 



accuracy, and establish the relationship of the effects found to familiar con- 
cepts in conventional nuclear physics, makes these degrees of freedom rather 
appropriate to look at the EMC and related effects. 
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Figures captions. 



• Figure 1: Selfenergy diagrams of the nucleon. 

• Figure 2: Electromagnetic form factors for the cases 
a) free nucleon, b) Fermi sea with B baryons. 

• Figure 3: (a) Feynman diagram for deep inelastic electron-nucleon scat- 
tering and (b), electron selfenergy diagram associated. 

• Figure 4: Electron selfenergy diagram accounting for electron-pion deep 
inelastic scattering. 

• Figure 5: Diagrams of the electron selfenergy including lph, lAh, lphlAh, 
etc.. 

• Figure 6: Two diagrams (a) direct and (b) crossed, which contribute to 
Compton 77r scattering. 

• Figure 7: Results of Rn(x) for 56 Fe. Solid line: using the spectral func- 
tion; dashed line: using the uncorrelated Fermi sea. Experimental points from 
ref. |32| (solids dots), ref. ||| (open squares). 

• Figure 8: Results of R(x) for 6 Li. Solid lines: whole calculation includ- 
ing the nucleons and the mesons; dashed line: contribution of the nucleons. 
Experimental points from ref. |M] (solid dots). Density for 6 Li from ref. [28 . 

• Figure 9: Same as fig. 8 for 12 C. Experimental points from ref. |M 
(solid dots), ref. |62| (open squares). 



• Figure 10: Same as fig. 8 for 40 Ca. Experimental points from ref. |)5 
(solid dots), ref. |3| (open squares). 

• Figure 11: Results for R(x) for 56 Fe. Solid line: whole calculation includ- 
ing the nucleons and the mesons; dashed line: contribution of the nucleons; 
dot-dashed line: contribution of nucleons plus pions. Experimental points 
from ref. [B^] (solid dots), ref. [p3| (open squares). 
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Table I 



<T> [MeV] < |E - M| > [MeV] \e A \ [MeV] 



6 Li 


UFS 


9.3 




22 


1 




5.2 




or 


lo.o 




oo 


Q 
O 




o.z 




UFS 


13.7 


(17.0) 


31 





(23.0) 


8.0 




SF 


31.7 


(37.0) 


50 


4 


(49.0) 


7.8 


m Ca 


UFS 


16.0 


(16.5) 


33 


5 


(26.6) 


8.5 




SF 


40.4 


(36.0) 


58 


6 


(52.1) 


8.6 




UFS 


16.1 


(17.0) 


34 


1 


(25.0) 


8.9 




SF 


40.4 


(33.0) 


58 


6 


(49.8) 


8.7 



UFS: uncorrelated Fermi sea. SF: spectral function. The numbers in brack- 
ets correspond to those obtained in ref. |54| respectively for the Hartree-Fock 



and spectral functions (the latter called there SRC, from short range correla- 
tions). 
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